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Abstract 

We extend chiral perturbation theory to study a meson gas out of thermal equi- 
librium. Assuming that the system is initially in equilibrium at Tj < T c and working 
within the Schwinger-Keldysh contour technique, we define consistently the time- 
dependent temporal and spatial pion decay functions, the counterparts of the pion 
decay constants, and calculate them to next to leading order. The link with curved 
space-time QFT allows to establish nonequilibrium renormalisation. The short-time 
behaviour and the applicability of our model to a heavy-ion collision plasma are also 
discussed in this work. 
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1 Introduction 



The chiral phase transition plays a fundamental role in the description of the plasma 
formed after a relativistic heavy- ion collision (RHIC), where it is imperative to use meson 
effective models to describe QCD. Two of the most successful approaches are the 0(4) 
linear sigma model (LSM), valid only for Nf = 2 light flavours, and Chiral Perturbation 
Theory (ChPT), based on derivative expansions compatible with the QCD symmetries, 
and whose lowest order action is the nonlinear sigma model (NLSM) [|TJ. In ChPT, the 
perturbative parameter is p/A x , with p a meson energy (like masses, external momenta 
or temperature) and A x ~ 1 GeV. Every meson loop is 0(p 2 /A^) and all the infinities 
coming from them can be absorbed in the coefficients of higher order lagrangians [1], 01. 

In thermal equilibrium at finite temperature T, the chiral symmetry is believed to be 
restored at T c ~ 150-200 MeV ||. In fact, near T c , the mean-field LSM is well known to 
undergo a second-order phase transition. The NLSM is equally valid for reproducing the 
phase transition, provided one works in the large- iV limit |4j]. Strictly within ChPT, the 
low-temperature meson gas has been studied |5j, |6| on expansion in T 2 /A 2 , predicting 
the correct behaviour of the observables as T approaches T c . 

The equilibrium assumption is not realistic if one is interested in the dynamics of the 
expanding plasma formed after a RHIC, where several nonequilibrium effects could be 
important. One of them is the formation of disoriented chiral condensates (DCC), regions 
in which the chiral field is correlated and has nonzero components in the pion direction 
JFj. As the plasma expands, long- wavelength pion modes — propagating as if they had 
an effective negative mass squared — can develop instabilities growing fast as the field 
relaxes to the ground state, an observable consequence being coherent pion emission ||. 
This issue has been extensively studied in the literature, mostly within the LSM assuming 
initial thermal equilibrium at Tj > T c , either encoding the cooling mechanism in the time 
dependence of the lagrangian parameters [8-11] or describing the plasma expansion in 
proper time and rapidity |T2j . This phenomenon has also been studied using Gross-Neveu 
models fl~3| . Another important nonequilibrium observable is the photon and dilepton 



production fll4 |, to which the anomalous meson sector could significantly contribute |10 



In this work we will construct an effective ChPT-based model to describe a meson gas 
out of thermal equilibrium, as an alternative to the LSM approach. Our only degrees of 
freedom will be then the Nambu-Goldstone bosons (NGB) and we will consider the most 
general low-energy lagrangian compatible with the QCD symmetries. We will restrict here 
to Nf = 2 (where the NGB are just the pions) and to the chiral limit (massless quarks), 
which is the simplest approximation allowing to build the model in terms of exact chiral 
symmetry. One of the novelties of our approach is to exploit the analogy between ChPT 
and the physical regime where the system is not far from equilibrium and then a derivative 
expansion is consistent. 



2 The NLSM and ChPT out of equilibrium 

We will take the system in thermal equilibrium for t < at a temperature T; < T c and 
for t > we let the lagrangian parameters be time-dependent. We are also assuming 
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Figure 1: The contour C in complex time t. The lines G\ and C2 run between t{ + it 
and tf + it and tf — it and tj — ie respectively, with e — > + . 

that the system is homogeneous and isotropic. The generating functional of the theory 
can then be formulated in the path integral formalism, by letting the time integrals run 
over the Schwinger-Keldysh contour C displayed in FigJI] [15-18]. We will eventually let 
U — > —00 and tf — > +00, although we will show that our results are independent of U and 
tf. We remark that, even in that limit, the imaginary-time leg of C has to be kept, since it 
encodes the KMS equilibrium boundary conditions [I7|, [18], [19] . With these assumptions, 
our low-energy model will be the following nonequilibrium NLSM 

S[U] = J c d\ tr djj\x, t)d»U{x, t) (1) 

where J c d 4 x = J c dt J d 3 x, U (x, t) G SU(2) is the NGB field, satisfying U (x, U + ifii) = 
U{x,U) with fa = Tr 1 , and f(t) is a real function which in equilibrium and to the 
lowest order (see section |j) would be / = f n ~ 93 MeV (the pion decay constant) i.e, 
f(t<0) — f. Note that f(t) cannot be analytic at t — and, in particular, it could 
be discontinuous, like the meson mass in quenched LSM approaches [8-11]. This is a 
consequence of the nature of our approach, since the system is driven off equilibrium 
instantaneously. An alternative, which we will not attempt here, is to choose f(t) analytic 
Vt, having equilibrium only at t = U [p0[| . Thus, the temporal evolution of our results 
will start at t = + , an infinitesimally small response time. As we will see below, our 
approach is consistent because the discontinuities at t = appear to NLO. 
We will parametrise the field U as 

U(x, t ) = jrr ) { [f^) - t)} V2 1 + ^vr a (f , t)} ; a = 1, 2, 3, (2) 

where 7r 2 = n a n a , ir a the pion fields satisfying ir a (ti + ifij) = 7r a (ti) and / and r a are 
the identity and Pauli matrices. Note that with the choice (§) we recover the canonical 
kinetic term in the action after expanding U in powers of n. Other choices amount to 
a time-dependent normalisation of the pion fields and should not have any effect on the 
physics (see section For instance, if we redefine 7r a = n a f '(0) / / (t) , the action for the 
7r a (x, t) fields is just the equilibrium NLSM multiplied by the time-dependent scale factor 
f (t)/f 2 (0) (see below). 
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Our action ([I]) is manifestly chiral invariant (U(x) —> LU(x)R)). Notice that we work 
in the chiral limit and hence there are no explicit symmetry-breaking pion mass terms 
in the action. The conserved axial and vector currents for the chiral symmetry can be 
derived by applying the standard procedure |E 0], so that the axial current reads 

A«(x,t) = itP-tv\r a (u^d^U -Ud,U^} (3) 



Let us now discuss how to establish a consistent nonequilibrium ChPT. The new 
ingredient we need is the temporal variation of f{t). We will then consider 
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(4) 



and so on, the rest of the chiral power counting being the same as in equilibrium. There- 
fore, in our approach we treat the deviations of the system from equilibrium perturba- 
tively, following the ChPT guidelines. Thus, we will expand our action (U) to the relevant 
order in pion fields and take into account all the contributing Feynman diagrams. The 
loop divergences should be such that they can be absorbed in the coefficients of higher 
order lagrangians, which in general will require the introduction of new time-dependent 
counterterms (see below). Notice also that according to (f|), we can always describe the 
short-time nonequilibrium regime, just by expanding f(t) around t = + . In fact, for times 
t < f~ x , that is equivalent to a chiral expansion, since then /(0 + )i//(0 + ) = 0(p/A x ) and 
so on. Nonetheless, we stress that the conditions (§) do not imply working at short-times, 
but just to remain close enough to equilibrium. 

To leading order in 7r fields, the action (P, after using @, reads 
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2 Jc 



fit) 



S [tt] = -- / (fxTi a (x,t) D + m 2 (t) 7r a (x,t) with m 2 (t) = -77^, (5) 



/(*)' 



where we have partial integrated in C. Thus, the leading order nonequilibrium effect of 
our model can be written as a time-dependent pion mass term, which, as commented 
before, is a common feature of nonequilibrium models [9-12]. Notice that m 2 (t) can be 
negative, so that our model accommodates unstable pion modes, whose importance we 
have discussed before. Note also that this mass term does not break the chiral symmetry, 
i.e, the axial current is classically conserved. Indeed, to leading order we have, from @, 

A a ,ix, t)} L ° = -f(t)d,7c a (x, t) + S, fit)n a (x, t) , (6) 

which satisfies d^A, = using [□ + m 2 (t)] ix a = 0, the equations of motion to the same 
order. Had we included the pion mass term — explicitly breaking the symmetry — the 
instabilities threshold, to leading order, would have been m 2 (t) < —m 2 instead. 

It is very interesting to rephrase our model as a NLSM in a curved space-time back- 
ground g, u , which reads [2! 



S 9 [U] = ^- J c d 4 x( v ^)g^trd,U\x)d u Uix)+£S R [U,R} (7) 
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plus U independent terms, where g = det g and the last term accounts for possible cou- 
plings between the pion fields and the scalar curvature R(x) (like R{x)<p 2 for a free scalar 
field pl|). Now, notice that our nonequilibrium model ([l|) is obtained by writing U(x) 
in the tt parametrisation discussed before (i.e, with f(t) replaced by /(0) in (^)), choos- 
ing £ = (minimal coupling) and a spatially flat Robertson- Walker (RW) space-time 
in conformal time, whose line element is ds 2 = a 2 (r])[dr] 2 — dx 2 }, with the scale factor 
a (v) = f{v)/f(ty- Our effective theory is then not only suitable for a RHIC environment, 
but also in a cosmo logical framework. Notice also that if we take £ 7^ 0, the lowest order 
Sr term we can construct has the form of an effective mass term breaking explicitly the 
chiral symmetry. In fact, it is not difficult to see that we could cancel the m 2 (t) term in ([5]) 
by choosing £ = 1/6, which is the value rendering the theory scale invariant pi| . This is 
just a consequence of the lagrangian chiral and conformal symmetries being incompatible 
in a curved background |22] or, equivalently, at nonequilibrium. In other words, for £ = 

we may 



-which is our choice, since we want to preserve chiral symmetry, as in 22 



interpret the m 2 (t) term, in the chiral limit, as the minimal coupling with the background 
yielding chiral invariance. 

The above equivalence turns out to be very useful to renormalise our model, consis- 
tently with ChPT. In fact, all the one-loop divergences arising from (0) can be absorbed 
in the coefficients of the 0(p A ) action 5*4, which consists of the Minkowski terms with 
indices raised and lowered with g^ plus new chiral-invariant couplings of pion fields with 
the curvature 



22 . In the chiral limit, those new terms read 



3? M 



d 4 x 



c 



g) [L lx R{x)g» v + L 12 R^} trd ll U\x)d u U{x) 



- / d A X7T a 
2 Jc 



fi(t)d 2 - f 2 {t)V 2 + m\{t) 



7T 



(8) 



where R^ v is the Ricci tensor, L\\ and Lyi are two new low-energy constants |22] and we 
have given the two-pion contributions in the parametrisation (g), after partial integration, 
with our RW metric, where 



h(t) 
h{t) 

m\(t) 



12 



(6Ln + L 



12; 



; Pit) 
f(t) 



+ L 



fi(t)f(t) + fi(t)f(t) 
f(t) 



12' 



[f(t)r 



(9) 



for t > and fi(t < 0) = 0. The above terms are the only ones in S4 containing two pions 
and they will renormalise purely nonequilibrium infinities — which are time-dependent 
and vanish for t < — . It is important to bear in mind that to cancel the one-loop new 
divergences only Ln needs to be renormalised, whereas L12 
back to this point below. 



-^12 



22|. We will come 



Next, we will concentrate on the Green functions time-ordered along C WA ITS] . Unless 



otherwise stated, we will be using the parametrisation (0) in the remaining of this work. 
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The two-point function defines the pion propagator G ab (x,y) = — i < TcTi a (x)TT a (y) >, 
which to leading order G^i^x, y) = S ab Go(x,y), by isospin invariance, and 



[n x + m 2 (x°)}G Q (x,y) = -8 c (x° - y°)5®(x - ft 



(10) 



with KMS equilibrium conditions G^(x,ti — if3i',y) = GQ-(x,ti;y), the advanced and 
retarded propagators being defined as customarily along C. Notice that G(x,x') = 
G(t,t',x — x') due to the nonequilibrium lack of time translation invariance. There- 
fore, we will define, as customarily, the "fast" temporal variable t — t' and the "slow" one 
t = (t + f)/2, so that F(q ,u q ,r) and F(u> g ,t,t'), with u 2 = \q\ 2 , will denote, respec- 
tively, the fast and mixed (in which only the spatial coordinates are transformed) Fourier 
transforms of F(x, x'). Note that F(q , w q , r) depends separately on q and u q because of 
the thermal loss of Lorentz covariance and has the extra nonequilibrium r-dependence. 
Then, in the mixed representation, fllOD becomes 



^+-;+- 2 (t) 



G Q (u q ,t,t') = -5 c (t-t') 



The general solution of (|TT1) is only known explicitly for some particular choices of 
m 2 (t) pl| , p~T| , 0. Formally, we can write it as a Schwinger-Dyson equation as 



G Q (aj q , t, t') = G e q (uj q , t-t')+ / dzm 2 (z)G e q (uj q , t - z)G (uJ q , z, t') 



c 



(12) 



with GQ q (u> q ,t — t') the equilibrium solution of (|TT|) , i.e, with m 2 {t) = 0. 

Another object of interest for our purposes is the Lehman spectral function p(x, y) = 
G > (x,y)—G < (x,y) |16|, which in equilibrium to leading order is p eq {q) = — 27risgn(g )<5(<? 2 ) 
[ |r§| . Note that, by construction, G > (x,y) = G < (y,x), so that p(x,y) = —p(y,x) and 
p(qo,u q ,T) = —p(—q ,uj q ,r). The normalisation of po is 



1 

2vri 



+oo 



q Po(q ,Ug,T) 



dpo(u q ,t,tf) 



dt 



-1, 



(13) 



t=v 



which can be readily checked by using (1T3) and p (u q ,t,t) = 0. 



3 Next to leading order propagator 

We will now obtain the NLO correction to the propagator. For that purpose, we need the 
action in ([D up to four-pion terms: 



S[*] 



SoM + -j c d*xt [W) 



+ - vr 



2\2 



Id) fHn 
J(t) f 2 (t), 



d^ a d»n b ir a ir b 
+ 0(7r 6 )) 



(14) 



plus the two-pion terms in (|8|). The two diagrams contributing are, respectively, a) and 
b) in Fig.||. 
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b) 




Figure 2: Diagrams contributing to the NLO pion propagator (a,b) and axial-axial 
correlator (c). The black dot in b) represents the interaction coming from Sj± in (H) 



Let us concentrate on Gn (i.e, t, t' G C\ in Fig.|l|) and t and t' positive, which is the 
relevant case for our purposes, as commented above. We will use dimensional regularisa- 
tion (DR), so that evaluating the above diagrams, using (0) with (JW(0) = 0, and after 
some algebra, we obtain in the mixed representation to NLO 



G>(t,t') = G> n (t,t')[l--[f 1 (t) + f 1 (t l )}) + 



2Un 



coth 



12/- 



COS [cJ q (t + t')] 



+ i 



dt 



dt 



dt [A^t, cu q )G>(t, t)G>(t, t') + A 2 (t)G>{t, t)G>(t, t' 
Ax(f, u q )G<(t, t)G<(t, f) + A 2 (t)G<(t, t)G<(t, t' 
Ai(t, u q )G<(t, t)G>(t, + A 2 (t)G<(t, t)G>(t, f 



(15) 



and Gf x {t : t') = G xl (t',t), where we have suppressed for simplicity the u q dependence of 
the propagators, the dot denotes d/dt, 



Ai(i,u q 



Pit) 



- ^ Go® " 2G (t) + 4^|g (*) 



+ ^ 



h{t)-h{t) 



— % 



mm , i f - 

~75T + 2 /l(t) 



A 2 (t) 



Gog) 
/ 2 (*~) 



(16) 
(17) 



and Gq(z°) = Gq(z, z) is the equal-time correlation function. We observe that (|T5|) is ti 
and independent, which is a good consistency check. Notice also that by replacing the 
equilibrium propagators in (|15D , we recover 



Git(t - 1') = gsS (t-o i 



T 2 



12/2 



which agrees with [|j] (note that we have derived it for the contour C, including both 
imaginary-time and real-time thermal field theory) and is finite in the chiral limit, where 
there is no tadpole renormalisation in DR 0. However, out of equilibrium, the NLO 
propagator is in general divergent, even in the chiral limit, and the infinities have to be 
absorbed in the two-pion counterterms in (pi). 
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4 The nonequilibrium pion decay functions 



In a thermal bath, the concepts of LSZ and asymptotic states are subtle, and so is then the 
extension of low-energy theorems like PCAC. Thus, pion decay constants are more conve- 
niently defined through the thermal axial-axial correlator A°? u (x, y) =< TcA a ^(x)A b ^(y) >. 
At T ^ the loss of Lorentz covariance in the tensorial structure of A^ v implies that one 
can define two independent and complex (spatial) and /* (temporal), their real and 
imaginary parts being related respectively with the pion velocity and damping rate in the 
thermal bath [p3|| . Nevertheless, to one-loop in the chiral limit one has || f|, ^3] 

[f:(T)} 2 =[fl(T)] 2 = / 2 (l-^, (19) 

with T c = y/Ef^ — 228 MeV. Despite it being just the lowest order in the low temperature 
expansion, (|i~9"D predicts the right behaviour and a reasonable estimate for the critical 
temperature, although, strictly speaking, /^(T) is not the order parameter ||. To higher 

orders, and Im#V OP- 



Let us then analyze A ab u in our nonequilibrium model. The relevant quantity, as far as 
f n is concerned, is the spectral function p^ u = A> — A< u , with A°? v = 5 ab A^. We readily 
realise that p fMU ( ( lo,(l, T ) = ~Pu^(—qo, ~^ T )- Then, from rotational symmetry, 

Pij(q , q, t) = qiqjp L (qo, u q , r) + 5 ij p d (q , uj q , r) (20) 

with pL,d(<?o) = -pL,d{-qo) Q and p j0 (qo,q, T ) = qjPs(qo,u q ,r). Therefore, p^ is charac- 
terized, in principle by the four functions pi, Pd, ps and poo- However, they are related 
through the A^ conservation Ward Identity (WI) d^p fll/ (x,y) = d^p^ u (x,y) = 0, which 
also holds in our model. Thus, we get 

q°Poo(qo, UJ q , r) - UJ^p s (q , UJ q , r) - -p 00 (qo, u q , t) = 

q p s (q ,uj q ,T)-UJgp L {qo,uJ q ,T) + -ps{qo,uJ q ,r)+pd{qo,uJq,r) = 0, (21) 

where the dot denotes d/dr. Thus, only two components of p^ are independent, as in 
equilibrium 0, where there are no time derivatives in the above equation. 

At T = one has pl = 27r/^sgn(g°)<5(g 2 ), since there exist NGB states. That is not 
the case at T ^ 0, where the pion dispersion relation is not in general a 5-function Q . In 
fact, to define properly f n (T) requires taking the uo q — > + limit, in which a zero-energy 
excitation still exists [|j], although to NLO there is no need to take that limit. Extending 
these ideas to nonequilibrium, we will define the time-dependent pion decay functions 
(PDF) as 

2 1 f°° d 

[Kit)} = 7T lin \ / dq q p L (qo,u q ,t) = lim % — p L (u q , t, f) (22) 

Z7T LL) q ^0+ J-oo ui q ^0+ dt , 

1 f°° 

/*(*)/*(*) = — lim / dq p s {qo,u q ,t) = lim p s (u q , t, t) (23) 

% f°° 1 d 

fn(t)9n{t) = -— lim / dq q ps{q ,uj q ,t) = - lim — p s ((J q , t, t' 

llX oj„^0+ J-oo Z uj„^0+ at 



(24) 



3 Our and pd correspond in in the notation of to sgn(g°) p^Qo I ^qQ 2 and sga(q°)p^ respectively. 
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The functions f^(t) and Pit) are the nonequilibrium counterparts of the spatial and 
temporal pion decay constants respectively, whereas g w (t) vanishes in equilibrium. How- 
ever, the above PDF are related through the WI. Integrating in go in (0); we get 

mgAt) = l^lftwm], (25) 

so that only two PDF are independent, as in equilibrium [p3| . Let us now check the con- 
sistency of our definitions to leading order. From (|B|), p^iojq, t, t') = if(t)f(t')po(uj q , t, t') 
and pa = 0, so that, using (|13|) yields f^{t) 2 LO = f 2 (t), i.e, the PDF coincides with 
fit) to leading order, as it should be. Similarly, we find p{t)LO — f^if)hO — fit) and 
9i?(t)LO — f(t), so that our definitions are consistent to leading order. 
To NLO, we need the axial current up to 0(tt 3 ). From (||), 

A a »(x,t) = [A^t)] L ° - -L_ (V«V 2 - n%n a - 6 ^j§j^ a ) + ( 26 ) 

with in (||). Thus, according to our chiral power counting, we have three types of 
NLO corrections to A^. The first is the NLO correction to the pion propagator we have 
evaluated in section BL coming from the product of the 0(n) terms above. The second is 
the product of the 0{n) with the C(vr 3 ), represented by diagram c) in Fig.||, and the third 
comes from the modification in A^ due to the action (|8|), which amounts to prefactors 
[1 + fi{t)] and [1 + f2{t)] in Aq and Aj respectively. Then, evaluating p^ to NLO, after 
using ( p~5|) (we take, without loss of generality, both t,t f 6 C\ and positive ) and (f22"D-(p3|), 
we finally arrive to 

[Pit)} 2 = f 2 (t)[l + 2f 2 (t)-f l (t)}-2zG (t) (27) 
flit)] 2 = f(t)[l + f 2 (t)}-2zG (t) (28) 

for t > 0. This is the main result of this work. It provides the NLO relationship between 
the PDF and f(t) Q. Notice that Pit) ^ /*(t) to NLO, unlike the equilibrium case, due 
to the effect of nonequilibrium renormalisation. However, note that [/*(£)] 2 — [p{t)] 2 = 
f 2 {t) [f2(t)— fi(t)], which is finite, since it depends only on L 12 , which does not renormalise. 
This is indeed an interesting consistency check, because the one-loop infinities appearing 
in Go(t) can then be absorbed in Ln, rendering both /^(t) and Pit) finite. We also 
remark that we did not need to take u q — > + in (|22"|)-(|23D to arrive to (|27|)-(p8l) (there are 
still NGB to NLO). Note also that both are real to this order. We have performed the 
following consistency checks on (p7j)-(|28|): first, the equilibrium result ( |19D is recovered 
(for the contour C) simply by replacing Gq Q = — iT 2 / 12 and fi = f 2 = 0. Second, by 
calculating g n (t) from p s , through (|24|), we check explicitly that the WI (^5|) holds and, 
third, we have calculated A^ in the 7r parametrisation, arriving to the same result. 

Therefore, (|27D-(p8D allow to express nonequilibrium observables (like decay rates, 
masses, etc) to one loop in ChPT, in terms of the physical f n {t), which could be mea- 
sured, for instance, in nonequilibrium lepton decays tt — > lv\. At this stage one can follow 
different approaches. Exact knowledge of fit) would require to solve self-consistently the 

4 G (t), hit) and f 2 (t) depend implicitly on f(t), through (0) and (|). 
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plasma hydro dynamic equations or, equivalently, Einstein equations for the metric. Alter- 
natively, one can treat f(t) as external — so that Q27|)-(]28D provide the system response- 



and study simple choices consistently with (||) [pOfl . In what follows, we shall take fit) 
arbitrary and expand it near t = + , analysing thus the short-time evolution. 

For short times, the particular form of f(t) is not important and we can parametrise the 
nonequilibrium dynamics in terms of the values of f(t) and its derivatives at t = + . As we 
discussed in section 0, this approach is justified for times t < t max with t max ~ 1/ f n {0) ~ 
2 fm/c (compare to the typical plasma time scales 5-10 fm/c [12|]). The general solution of 
flTT|) with KMS conditions at ti can be constructed in terms of two independent solutions 
to the homogeneous equation, which have to be continuous and differentiable Vt G C so 
that the solution is uniquely defined ]I7|]. Therefore, they have to match the equilibrium 
solution and its first time derivative at t = 0. With these conditions and expanding both 
fit) and the solutions near t = + we find to the lowest order 



Gq (w„M) 



2u n 



coth 



1 



2.2 

m t 



0(mH 



(29) 



for t > 0, with m 2 = — /(0 + )//(0 + ). For m 2 < we see the unstable modes threshold, 
making the pion correlation function grow with time. The effect of those modes is not 
important for short times though, where the exponential growth of the correlator is not 
appreciable. Observe that in fl2"9| ) the time dependence factorises, so that the momentum 
dependence is the same as in equilibrium and then we can integrate it in DR, yielding the 
finite answer (ITS). Then, from (|2"7|)-(|2"g|) we get 



#*(*) 



rs,t 
JR 





0(P 3 /Al 



(30) 



for t > 0, where H = /(0 + )//(0), with the renormalised constants 



[f. 



Ri 



/ 2 (0+)+4 



' L 



12 



6L n W - L 12 H 2 



fk 



/ 2 (0 + )+4 -(Li 2 + 6L n )m 2 + L 12 H< 



(31) 



and where the H and m parameters (which are 0(p) and play the role of the Hubble 
constant and the deceleration parameter in the Universe expansion) also get renormalised, 
in terms of /i,2(0 + ), A,2(0 + ) and so on, but those are subleading contributions. Thus, 
for short times, all the effect of the S4 terms, which is Ti independent, is to redefine 
/?r(0 + ), since there are no infinities coming from Go in DR. We insist that this is just the 
effect of truncating the series in t and it is not true in general. Notice that ( |31]) implies 
necessarily a nonzero jump Af = f(0 + ) — / (see our comments in section 0) so that the 
divergent part of L\\ can be absorbed in / 2 (0 + ) rendering a finite A/jj* = A/^' 4 . In fact, 
that effect is very small compared to the other contributions in (|3"0"D, and so it is the 
difference f%(t) — fl(t), since L r n , V 12 — 10~ 3 p2| . Notice also that for the particular case 



H 



m 



> 0, renormalising as / 2 (0 + ) = f 2 + 24m 2 Ln we get = /* and A/^ =0. 
Finally, we will estimate some physical effects related to f n (t). For that purpose, we 
will ignore, for simplicity, the effect of L n and L 12 and, based upon (O), define the plasma 
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effective temperature as f 2 (t) = f 2 [l — T 2 (t)/T 2 ]. Therefore, we can also define a critical 
time as T(t c ) = T c and a freezing time T(tf) = 0. Thus, we will impose < T{t) < T c 
and then, through our short-time results for f n (t), determine either t c or tf, depending 
on the initial conditions (T(t) is just quadratic in time to this order). Notice that we are 
following a similar approach as in equilibrium when one extrapolates (O) until T = T c . 
Let us then take typical values Tj, \H\, \m\ ~ 100 MeV and retain only the leading order in 
x = Tf /T c 2 , consistently with the chiral expansion. Then, if H = 0, the system cools down 
until t 2 m 2 ~ —x(l + x) (t/ — 0.2 fm/c) if m 2 < 0, whereas for m 2 > it is heated until 
t 2 m 2 — 1 (t c ~ 2 fm/c), independent of Tj. For H > 0, there is cooling until £/|-ff| — x/2 
(t/ ~ 0.2 fm/c). Finally, for < and m 2 > there is heating until t c \H\ ~ (1 — 3x/4)/2 
(t c ~ 2 fm/c), whereas if m 2 < 0, there is heating until a maximum t m \H\ ~ (1 + x/2)/2 
and then cooling down until t/|-f^| — 1 + x (tf ~ 2.3 fm/c). We observe that the effect 
of the unstable modes (m 2 < 0) is always to cool down the system and that the freezing 
time for H < is much longer than that for H > 0. Some of these time scales are indeed 
longer than those to which our short-time approximation remains valid, but they have to 
be understood as estimates, similarly to estimating T c at equilibrium through ([HJ), even 
though the low T approach is less reliable near T ~ T c . 

Comparing with flT2"|| , naively identifying the LSM order parameter v(t) ~ fn{t) (in 
proper time), we see a similar short-time evolution, although our estimates for the time 
evolution duration are somewhat lower. This was expected, since the initial values in W2 



correspond to Tj ~ 200 MeV and \H\ ~ 400 MeV, which are too high for our low-energy 



approach. An important remark is that in typical simulations like |T2[], v{t) reaches a 
stationary value, about which it oscillates (thermalisation). It is clear that we cannot 
predict that type of behaviour only within our short-time approach, quadratic in time, 
but only estimate the time scales involved — similarly as to why ChPT cannot see the 
phase transition — . Therefore, in view of the above estimates, we believe that our ChPT 
model may be useful for studying the different nonequilibrium observables evolution, 
from a stage where some cooling has already taken place onwards. In principle, we could 
approach closer to T c by considering enough orders in our ChPT, although in practice, 
beyond one-loop, some resummation method, like large N, will need to be implemented. 



5 Conclusions and Outlook 

We have extended chiral lagrangians and ChPT out of thermal equilibrium. The chiral 
power counting requires all time derivatives to be 0(p) and to lowest order our model 
is a NLSM with / — > f(t). This model accommodates unstable pion modes and corre- 
sponds to a spatially flat RW metric in conformal time with scale factor a(t) = f{t)/f(0) 
and minimal coupling. We have exploited this analogy to establish the renormalisation 
procedure, which allows to construct the fourth order lagrangian absorbing all the loop 
divergences, which in general will be time-dependent. 

We have applied our model to study the time-dependent pion decay functions, ex- 
tending the equilibrium pion decay constants. In general there are two independent PDF, 
as in equilibrium, and to NLO in ChPT they already differ, unlike equilibrium, due to 
renormalisation. We have obtained them to NLO in terms of the equal-time correlation 
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function, analysing their lowest order short-time coefficients and their dependence with 
Tj, and discussing the relevant time scales involved within the context of a RHIC plasma. 

Among the aspects of our model which are worth to be studied further are the long-time 
evolution, by choosing suitable parametrisations for /(£), including the analytic approach, 
and the behaviour of the two-point correlation function at different space points, which 
would allow us to investigate the formation of regions of unstable vacua (DCC) f20|- Other 
applications and extensions, to be explored in the future include photon production in 
the pion sector (by gauging the theory and including 7r° anomalous decay), the quark 
condensate time dependence (by including the mass explicit symmetry-breaking terms), 
the Nf = 3 case, large N resummation and proper time evolution. 
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